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QUALITATIVE  MATRICES: 

STRONG  SIGN-SOLVABILITY  AND  WEAK  SATISFIABILITY 

VICTOR  KLEE  and  RICHARD  LADNER 

Introduction 

The  study  of  "qualitative"  linear  systems  was  suggested  by  Paul 
Samuelson  in  1947  in  his  influential  book  on  Foundations  of  Economic 
Analysis  [29].  Since  then  it  has  been  carried  on  by  economists,  mathe¬ 
maticians,  ecologists  and  chemists,  with  emphasis  on  questions  of  quali¬ 
tative  solvability  and  qualitative  stability.  Survey  articles  have  been 
written,  and  references  collected,  by  Maybee  and  Quirk  [24],  Maybee  [23] 
and  Quirk  [28].  The  present  paper  is  devoted  to  a  rather  general  form 
of  qualitative  solvability,  of  which  strong  sign-solvability  is  a  spe¬ 
cial  case,  and  to  a  closely  related  notion,  weak  satisfiability,  from 
propositional  logic. 

The  simplest  qualitative  conditions  are  those  involving  sign-patterns. 
For  a  (real)  m  x  n  matrix  A  ■  (a^ )  ,  let  Q(A)  denote  the  convex  cone 
consisting  of  all  matrices  that  have  the  same  sign-pattern  as  A  .  Thus 
a  matrix  B  ■  (b^j)  belongs  to  Q(A)  if  and  only  if  B  is  an  m  *  n 
matrix  with  sgn  b^  -  sgn  a„  for  all  1  and  j  .  Here  sgn  t  Is  -, 

0  or  +  according  to  whether  the  real  number  T  is  <  0,  ■  0,  or  >  0  . 

Points  of  R®  are  taken  as  column  matrices  unless  the  contrary  is 
specified.  For  c  e  R®  ,  (Arc)  will  denote  the  linear  system  whose  coef¬ 
ficient  matrix  la  A  and  "constants  colon"  Is  c  .  The  system  (Arc) 

Is  said  to  be  eign-solvable  If  A  Is  square,  the  system  Is  solvable,  and 
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both  its  solvability  and  the  sign-pattern  of  its  solution  depend  only  on 

the  sign-patterns  of  A  and  c  .  That  is,  there  exists  x  vith  Ax  *  c 

such  that  for  each  B  e  Q(A)  and  d<Q(c)  it  is  true  that  (i)  there 

exists  y  with  By  ■  d  and  (ii)  for  each  such  y  ,  Q(y)  -  Q(x)  .  The 

system  (A:c)  is  strongly  sign-solvable  if  there  exists  x  as  described 

with  all  coordinates  of  x  different  from  0  . 

As  is  shown  in  Section  1,  the  study  of  strong  sign-solvability  is 

easily  reduced  to  the  study  of  what  we  call  S-matrlces.  An  S-matrix  is 

nH"l 

an  m  *  (m  +  1)  matrix  A  such  that  for  each  B  e  Q(A)  and  x  e  R 

with  Bx  ■  0  ,  all  coordinates  of  x  are  of  the  same  sign;  equivalently, 

each  such  B  has  as  its  nullspace  a  line  that  Intersects  the  open  posi- 
flx4*l 

tlve  orthant  of  R  .In  1962  Lancaster  [18]  suggested  a  general  form 
for  S-matrices;  more  inclusive  forms  were  then  described  by  Gorman  [12] 
and  Lancaster  [19].  In  1965  Lancaster  [20]  initiated  the  algorithmic 
approach  to  S-matrices,  based  on  duality  properties  of  convex  cones.  Two 
central  problems  have  been  those  of  finding: 

(a)  a  small  collection  of  ’’standard  forms"  of  S-matrices  to  which 
all  such  matrices  can  be  reduced  by  certain  elementary  trans¬ 
formations; 

(b)  a  fast  algorithm  for  the  recognition  of  S-matrlces. 

(The  algorithms  of  Lancaster  [20,  21]  are  of  exponential  cooplexity.) 

It  is  not  clear  to  us  that  one  should  expect  e  truly  useful  solution 
of  (a).  In  any  case.  Section  2  describes  some  ways  of  constructing  S- 
mat rices  and  the  closely  related  NW -mat rices.  A  conjecture  of  Corman  [12] 
Is  elucidated  but  not  settled. 
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In  Section  4  the  goal  (b)  is  attained  by  reducing  the  S-matrix  recog- 

2 

nit ion  problem  to  a  sequence  of  m  +  2m  +  2  linear  feasibility  tests. 

For  the  recognition  of  S-matrices,  per  se,  it  is  probably  more  efficient 
to  use  the  graph-theoretic  characterisation  of  sign-solvability  established 
by  Bassett,  Maybee  and  Quirk  [2]  and  discussed  further  by  Maybee  [22,  23]. 
However,  our  method  goes  considerably  beyond  the  case  of  S-matrlces,  apply¬ 
ing  whenever  the  columns  of  the  m  x  (m+  1)  matrix  A  are  restricted 
only  by  membership  in  given  polyhedral  cones.  If  the  feasibility  tests  are 
made  by  means  of  the  recent  Shor-Khachlan  method  [15,  31],  the  worst-case 
complexity  of  our  algorithm  is  bounded  by  a  polynomial  in  the  length  of  the 
binary  encoding  of  the  description  of  the  cones.  Though  the  Shor-Khachian 
procedure  is  not  well-suited  to  actual  computation  (see  the  comnents  in  Sec¬ 
tion  4),  this  result  is  of  theoretical  Interest  because  the  S-matrlx  recog¬ 
nition  problem  is  in  a  sense  close  to  the  boundary  of  NP-completeness . 

In  the  language  of  propositional  logic,  recognizing  that  a  given 
m  x  (m  +  1)  matrix  is  not  an  S-matrlx  amounts  to  recognizing  that  an  asso¬ 
ciated  Boolean  formula,  consisting  of  the  conjunction  of  n  +  1  disjunctive 
clauses  in  n  propositional  variables,  is  weakly  satlsfiable.  Although 

this  problem  can  be  solved  in  polynomial  time,  we  show  in  Section  3  that  for 

1/k 

each  k  >  0  the  problem  of  recognizing  the  weak  satisfiability  of  n  +  Ip  J 
clauses  in  n  variables  in  HP-complete. 

The  section  headings  are  as  follows:  §1.  Strong  sign-solvability, 
S-matrlces  and  weak  satisfiability;  i2.  Constructions  of  NW-matrices; 

13.  The  recognition  of  weak  satisfiability;  S4.  The  recognition  of  S- 
matrlces;  15.  Open  problems. 
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1.  Strong  Sign-solvability.  S-matrices  and  Weak  Satisfiability 


The  following  remark  Is  straightforward. 

THEOREM  1  If  c  is  a  column  of  an  S -matrix  A  and  B  is  the  square  matrix 
formed  by  A 's  other  columns  then  the  system  (B:c)  is  strongly  sign-solvable. 
Now  suppose,  conversely,  that  B  is  an  m  *  m  matrix,  the  system  (B:c)  is 
strongly  sign- solvable,  and  Bx  ■  c  .  Let  A  be  the  m  x  (m  +  1)  matrix 
whose  (m  +  1 )th  column  is  -c  and  whose  j th  column  is  ( for  1  <  j  <  m) 
the  jth  column  of  B  or  the  negative  of  that  column  according  as  x ^  >  0 
or  Xj  <  0  .  Then  A  is  an  S -matrix. 

In  view  of  condition  (11)  of  the  next  result,  we  may  regard  S  as  standing 
for  "simplex"  as  well  as  for  "strongly  sign-solvable".  We  could  avoid  (11), 
and  establish  the  equivalence  of  (1)  and  (111)  directly,  by  appealing  to  stand¬ 
ard  results  on  convex  polyhedral  cones  [9,  11).  In  fact,  that  was  done  by 
Lancaster  (20].  However,  we  want  to  Include  (11)  because  It  provides  useful 
geometric  Insight.  A  more  geometrical  formulation  of  (ill)  is  as  follows: 
for  each  open  halfspace  H  In  R™  whose  bounding  hyperplane  passes  through 
the  origin  there  Is  a  column  v^  of  A  such  that  the  cone  Q(Vj)  con¬ 
tained  in  H  . 


THEOREM  2  Suppose  that  A  is  an  m  *  (m  +  1)  matrix  and 

are  the  colwms  of  A  .  Then  the  following  three  conditions  are  equivalent: 

(1)  A  is  an  S -matrix; 

(11)  for  each  choice  of  Cj  e  Q(v1),...,cofl  e  <^vm+l*’  the  cj'e 
are  the  vertices  of  an  m-simplex  whose  interior  contains  the 
origin; 

(ill)  for  eaoh  nonzero  t  e  8®  there  exists  J  such  that  cl t  >  0 
for  all  c  «  Q(Vj) . 
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Proof.  Suppose  first  that  (i)  holds  and  let  e  Q(Vj)  for 
1  i  j  £  m  +  1  .  Let  B  be  the  m  *  (m  +  1)  matrix  whose  jth  column 
Is  Cj  .  Since  the  rank  of  B  Is  at  most  m  ,  there  exists  a  nonzero 
x  such  that  Bx  ■  0  .  By  (1),  each  such  x  has  exclusively  negative 
or  exclusively  positive  coordinates,  whence  It  follows  readily  that  the 

set  (c^ . c[|t>  Is  affinely  Independent  and  the  origin  is  a  strictly 

positive  convex  combination  of  the  c^'s  .  But  that  is  the  content  of 
(11)  .  Hence  (1)  Implies  (11) ,  and  the  argument  Is  easily  reversed  to 
show  (11)  Implies  (1) . 

If  (ill)  fails  there  exists  z  4  0  and  there  exist  Cj  e  Qfv^ . 

cnt)  ^  e  Qfv^  ^)  such  that  c^z  £  0  for  all  j  .  Plainly  that  contradicts 
(11),  and  hence  (11)  Implies  (ill).  Suppose,  finally,  that  (ill)  holds 
and  Cj  €  Q(Vj)  *or  1  5  J  Sm  +  1  .  Let  K  denote  the  convex  hull  of 
the  cj's  *  follows  from  (111)  that  K  Intersects  every  open  halfspace 

whose  bounding  hyperplane  passes  through  the  origin.  But  then  the  origin 
must  be  Interior  to  K  ,  and  since  K  has  at  most  m  +  1  vertices  It  must 
In  fact  be  an  m-slmplex.  That  shows  (ill)  implies  (11)  and  completes  the 
proof  of  Theorem  1.  □ 

Now  let  us  define  a  W -matrix  as  an  m  x  n  matrix  A  for  which  there 
exists  B  e  Q(A)  and  nonzero  x  e  Rn  with  Bx  i  0  (that  is,  each  coor¬ 
dinate  of  Bx  is  i  0) .  An  NW -matrix  is  one  that  is  not  a  W-matrix.  Con¬ 
dition  (111)  above  asserts  the  transpose  A*  is  an  NW-matrix.  In  terms 
of  the  notions  defined  in  the  next  paragraph,  W  may  be  regarded  as  stand¬ 
ing  for  "weakly  satlsfiable".  The  connection  with  satisfiability  in  propo¬ 
sitional  logic  is  clarified  in  Section  3. 


When  x  ■  (x^,...,xn)  and  y  -  (yp...,yn)  are  ordered  n-tuples 
of  real  numbers,  we  say  that  x  hits  <resp.  satisfies >  y  if  there 

exists  j  such  that  x^y^  4  0  <resp.  >  0> ;  and  x  misses  y  if  x 
does  not  hit  y  .  Further  x  weakly  satisfies  y  if  x  satisfies  y 
or  x  is  nonzero  but  misses  y  .  And  x  weakly  satisfies  the  m  x  n 
matrix  A  if  x  weakly  satisfies  each  row  of  A  .  Note  that  A  is 
weakly  satisfiable  if  and  only  if  A  is  a  W-matrix. 
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2.  Construction  of  NW-matrices 


In  constructing  m  x  n  NW-matrices,  we  are  especially  interested  in 
those  for  which  m  ■  n  +  1  ,  since  they  are  the  transposes  of  S-matrlces. 

When  m  £  n  there  are  no  m  x  n  NW-matrices  because  each  set  of  n  points 
in  Rn  lies  in  a  closed  halfspace  whose  bounding  hyperplane  passes  through 
the  origin. 

If  A  is  an  NW-matrlx  then  so  is  every  matrix  obtained  from  A  by 
permuting  rows,  permuting  columns,  replacing  columns  by  their  negatives, 
and  changing  the  magnitudes  (but  not  the  signs)  of  individual  entries. 

These  operations  provide  the  natural  equivalence  relation  for  NW-matrices. 
Since  only  sign-patterns  are  Involved,  each  equivalence  class  may  be  repre¬ 
sented  by  various  arrays  of  the  symbols  -,  0  and  +.  We  take  as  the  canonical 
representative  the  one  which  is  lexicographically  first  when  each  matrix  is 
considered  as  the  sequence  formed  by  writing  down  its  successive  rows.  The 
lexicographic  ordering  is  based  on  the  ordering  -  <  0  <  +  of  the  sign- 
symbols.  Thus,  for  example,  —  and  +  are  the  two  2xl  NW-matrices  and 
+  is  the  canonical  representative  because  - +•  precedes  +  -  in  the  lexi¬ 
cographic  ordering.  Fig.  1  shows  the  two  canonical  representatives  of  the 


Fig.  1:  The  two  canonical  3x2  NW-matrlces  and  associated  2-simplices 
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By  means  of  a  computer  search  we  have  verified  that  there  are  precisely 
ten  equivalence  classes  of  A  x  3  NW-matrices.  Their  canonical  representa¬ 
tives  are  shown  in  Fig.  2. 

--+  --+  -0+  -0+  -  0  + 

-  +  0  0  +  0  -  +  0  0  +  -  0  +  0 

+  00  +00  +00  +00  +00 

-  --0  0  0  -  -  0 

00+  -0-  -  +  0  —  +  0  0  0- 

0+0  0  +  +  00-  +0-  0  +  0 

+  00  +00  +0+  +0+  +0  + 

Fig.  2:  The  ten  canonical  A  x  3  NW-matrices 

Each  of  the  3  x  2  classes  can  be  constructed  from  the  2x1  class 
in  the  manner  shown  in  Fig.  3,  and  all  but  the  next-to-last  A  x  3  class 
arises  from  a  3x2  class  in  the  same  way.  In  Fig.  3,  the  new  row  has 
all  entries  0  except  for  a  single  nonzero  s(-  or  +)  in  the  new  column 
has  s  in  the  new  row,  *",s  in  some  other  position,  and  its  remaining 
entries  are  all  0  or  “ '6  .  (We  agree  that  “1-  is  +  and  — '+  is  -  . ) 

n 

1 

Fig.  3:  Forming  an  (m  +  1)  x  (n  +  1)  NW-matrix 
from  an  m  x  n  NW-matrix  A 

With  m^  ■  m^  ■  k  ■  2  and  n^  ■  n^  ■  ■  1  ,  the  construction  shown 

In  Fig.  A  yields  a  A  x  3  NW-matrix  that  is  equivalent  to  the  next-to-last 
one  of  Fig.  2.  In  general,  one  starts  from  a  k  x  t  NW-matrix  C  and 
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Ca)  T  is  a  rooted  tree  in  which  each  internal  node  has  at  least 
two  sons; 

(b)  for  each  node  j  of  T,  Y(j)  is  a  nonempty  subset  of  X  ; 

(c)  for  the  root  r  of  T,  Y(r)  ■  X  ; 

(d)  if  i  is  an  internal  node  of  T  and  S(l)  is  the  set  of  all 

sons  of  1  in  T  then  (Y(j)  ;  ]  e  S(l)}  is  a  partition  of 
Y(i)  —  that  is,  Y(j)  n  Y(j')  ■  0  for  each  choice  of  distinct 
j,  j'  e  S(i),  and  ‘ ni>  • 

For  our  purposes,  no  generality  is  lost  by  adding  the  requirement  that 
X  is  a  finite  linearly  ordered  set  and  for  each  node  j  the  set  Y(j)  is 
an  interval  with  respect  to  the  ordering.  Fig.  5  depicts  a  tree  with  root 
r  *  1  and  node-set  {1,...,8},  Fig.  6  an  associated  partition-tree  for  the 
set  X  ■  {1, ... ,11}  . 


Fig.  5:  A  tree  with  root  1  and  node-set  {1,...,8} 
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Y(l)  Y(2) 

Fig.  6:  A  partition-tree  for 

associated  with  the  tree  of  Fig.  5 


THEOREM  3  Suppose  that  (T,  Y)  is  a  partition-tree  for  the  set  , 

I  is  the  set  of  all  internal  nodes  of  T  ,  and  L  is  the  set  of  all  leaves 
l  of  T  for  which  |Y(£)|  2  2  .  For  each  lei  let  be  an  |S(i) I  *  n± 
m-matrix  and  for  each  l  e  L  let  be  a  |Y(£)  |  *  NW -matrix.  Then 

the  following  construction  yields  an  m  *  (Ej£lULB.j  )  NW  -matrix  B 

(i)  For  each  l  e  L  let  B  have  colxtms  which,  in  the  rows 

corresponding  to  ~  Y (l)  < resp .  Y(£)>  have  all  entries 

0  (peep,  a  copy  of  A^); 

(ii)  For  each  lei  let  Ri»****R|s(i)  |  the  sets  Y(j)  for 

j  e  S(l)  .  Let  B  have  columns  which  have  all  entries  0 
in  the  rows  corresponding  to  {1, . . .  ,m}  ~  Y(i)  and  which  satisfy 
the  following  condition  for  1  s  k  i  S(l)  :  each  row  corresponding 


-12- 


to  has  all  entries  0  or  is  a  copy  of  the  k th  row  of  A^  ,  and 
there  is  ( for  eaah  k  )  at  least  one  such  copy. 

Before  proving  Theorem  2,  let  us  illustrate  its  use  in  the  construction 
of  NW-matrices.  To  obtain  the  construction  of  Fig.  3  from  Theorem  2,  let  T 
and  Y  be  as  in  Fig.  7,  let  be  +  or  +  ,  and  let  A2  be  the  A  of 

Fig.  3.  There  is  no  A^  because  |Y(3)|  ■  1  . 


Y(  1)  Y(3) 

Fig.  7:  Obtaining  Fig.  3’s  construction 
as  a  special  case  of  Theorem  2 


To  obtain  the  construction  of  Fig.  4  from  Theorem  2,  let  T  consist  of  a 
root  1  directly  joined  to  nodes  2,...,k  -f  1  ,  let  mi  *or 

1  s  j  S  k  ,  and  let  Y(l)  -  {l,...^},  Y(2)  -  {lt...,o1> .  Y(k  +  1)  - 

{ok_j  +  l,...,ok)  .  Let  the  Aj,  A2,  ....  A^  of  Theorem  2  be  respectively 
the  C,  A^,  . ...  A^  of  Fig.  4. 

For  one  more  illustration,  let  (T,  Y)  be  as  in  Figs.  5-6  and  let 

Ai  -  A2  -  +,  A6  -  A?  -  ±.A3  -  A,.  -  - A8-0  +  . 

Because  of  the  freedom  ellowed  In  (11)  in  forming  the  colons  of  B  cor¬ 
responding  to  Internal  nodes  of  T  ,  the  number  of  11  *  10  NW-matrices 
that  can  be  constructed  in  the  manner  of  Theorem  2  with  the  specified  T, 
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Y  and  A^'a  la 

<29  -  1)  (22  -  1)  <25  -  1)  (23  -  1)  <22  -  1)  (23  -  1)  -  6,999,552. 


In  particular,  if  no  uae  la  made  of  the  freedom  allowed  in  (11)  (that  Is, 
If  each  row  correapondlng  to  la  a  copy  of  the  kth  row  of  A^) ,  the 
matrix  of  Fig.  8  la  obtained. 


0 

0 

0 

-  0  - 

0 
0 
0 
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0  + 
+  -  0 

+  +  0 

Fig.  8:  An  11  *  10 


0 

0 

0 

0 

0 

+ 

+ 

+  + 

0  0 

0  0 

0  0 


0  + 

0  + 

0 
0 
0 

0 

+  0 

0  0 

0  0 

0  0 

0  0 


+  0 

0  0 

0  + 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 


0 

+ 

0 

0 

0 

0 

0 

0 

0 

0 


NW -matrix  constructed  by  the  method  of  Theorem  2 


Proof  of  Theorem  2.  Suppose  there  exists  a  matrix  B  which  is  weakly 
eatl8flable  even  though  It  1 a  generated  by  a  partition-tree  (T,  Y)  and 
by  NW-matrices  A^  In  the  manner  of  Theorem  2.  Among  all  auch  B  ,  con- 
aider  one  for  which  the  nixnber  of  nodes  of  an  associated  T  Is  a  minimum. 
With  n  ■  Ij£iuLnj  *  the  manber  of  columns  of  B  ,  let  x  •  (x^,...,xn) 
be  an  n-tuple  that  weakly  satisfies  B  .  Let  H  denote  the  set  of  all 
nodes  h  of  T  auch  that  x  hits  aome  row  of  B  In  the  portion  of  B 
that  la  generated  by  the  node  h  and  the  NW -matrix  A^  .  Then  H  la 
nonempty,  for  all  nonsero  entries  of  B  are  generated  In  this  way  and  no 
column  of  B  la  aero.  Among  all  members  of  H  ,  choose  hQ  closest  to  the 


root  r  of  T  In  the  sense  thst  the  psth  from  r  to  1»q  in  T 
Intersects  H  only  st  hQ  .  Let  TQ  be  the  subtree  of  T  thst  Is 
rooted  st  h^  snd  let  Bq  be  the  submstrlx  of  B  thst  is  sssocisted 
with  Tq  .  Then  Bq  Is  weekly  sstisfied  by  the  relevsnt  portion  of  x  , 
snd  since  Bq  Is  slso  genersted  in  the  msnner  of  Theorem  2  it  follows 
from  the  silnimslity  of  T  thst  Tq  •  T  snd  hp  -  r  .  If  r  is  s  leaf 
of  T  then  B  is  the  NW-mstrix  snd  an  lmnediate  contradiction  en¬ 

sues. 

Now  suppose  that  the  root  r  is  an  internal  node  of  T  ,  snd  let 
x*  denote  the  restriction  of  x  to  the  column  associated  with  r  . 

Since  x*  i  0  snd  A ^  is  an  NW-matrix,  there  exist  r'  and  k  such 
thst  (in  the  notation  of  (11)  of  Theorem  2)  r'  e  S(r),  1  <  k  $  |S(r)|, 

R^  -  Y(r') ,  and  x*  hits  but  does  not  satisfy  the  kth  row  of  A^  . 

Let  T'  denote  the  subtree  of  T  that  is  rooted  at  r'  and  let  B' 
denote  the  submatrix  of  B  that  is  formed  from  the  rows  of  B  correspond¬ 
ing  to  R^  snd  the  columns  of  B  associated  with  nodes  of  T'  .  Let  x' 
denote  the  restriction  of  x  to  these  columns.  From  the  fact  that  x 
weakly  satisfies  B  while  x*  hits  but  does  not  satisfy  the  k**1  row 
of  Ar  it  can  be  deduced  thst  x'  weakly  satisfies  B'  .  But  B'  is 
formed  from  T'  snd  its  associated  Y(j)'s  and  A^'s  in  the  manner  of 
Theorem  2,  contradicting  the  minimality  of  T  snd  completing  the  proof.  □ 

Gorman's  1964  conjecture  [12]  on  the  construction  of  S~matrlces  was 
rephrased  in  [18,  23].  In  our  terms,  it  amounts  to  saying  that  if  JJ  is 
the  class  of  all  NW-matrlces  that  have  exactly  one  more  row  than  columns, 
then  each  member  of  g  is  equivalent  to  one  that  is  generated  in  the  manner 


of  Theorem  2  by  a  partition-tree  (T,  Y)  and  a  collection  of  NW-matrlces 
Aj  such  that  T  la  binary  (each  internal  node  has  exactly  two  eons)  and 
each  Aj  la  ♦  or  +  .  Ve  have  not  settled  this  conjecture.  More  gener¬ 
ally,  let  P  denote  any  one  of  the  following  three  construction  procedures: 
(a)  combining  the  method  of  Fig.  2  with  that  of  Fig.  3;  (b)  restricting  the 
method  of  Theorem  2  to  binary  trees;  (c)  using  the  method  of  Theorem  2  with¬ 
out  restriction.  Then  we  are  unable  to  say  whether,  for  an  arbitrary  member 
H  of  |j  that  has  more  than  two  rows,  a  member  equivalent  to  M  can  be 
generated  by  applying  P  to  members  of  JJ  that  are  smaller  than  M  . 
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3.  The  Recognition  of  U-matrlces 

For  «n  m  *  n  matrix  A  ,  weak  aatlsf lability  can  be  tested  in 
0(3nmn)  steps  by  generating  each  of  the  3n  -  1  possibilities  for  the 
sign-pattern  of  nonzero  y  e  Rn  and,  as  each  one  is  generated,  testing 
to  see  whether  it  weakly  satisfies  each  row  of  A  .  Similar  procedures 
were  suggested  by  Lancaster  [20,  21]  for  similar  purposes.  Although 
these  algorithms  are  finite  and  very  easy  to  program,  they  are  practical 
only  for  small  values  of  n  .  However,  we  show  below  that  for  general  m 
and  n  ,  and  even  for  the  cases  in  which  m  ■  n  +  Ln^J  for  an  arbitrary 
fixed  k  >  0  ,  the  problem  of  recognizing  weak  satisfiability  is  NF-complete 
and  hence  algorithms  requiring  an  exponential  number  of  steps  may  be  un¬ 
avoidable. 

It  will  be  convenient,  for  the  rest  of  this  section,  to  change  the 

language  from  that  of  matrix  theory  to  that  of  propositional  logic.  The 

n  columns  of  an  m  *  n  matrix  A  ■  (a^)  correspond  to  n  propositional 

variables  the  m  rows  of  A  correspond  to  clauses  Cp..., 

C  in  the  literals  ,  u u  ,  where  u.  stands  for  ~*u..  The 

m  i  n  i  n  J  j 

matrix  A  is  represented  by  a  Boolean  formula 

B  -  Cj  A  C2  A  ...  AC# 

in  conjunctive  normal  form,  where  the  clause  is  the  disjunction  of 

literals  obtained  from  the  i**1  row  of  A  as  follows: 

Uj  appears  in  if  and  only  if  a^  >  0  ;  u^  appears  in  if 
and  only  if  a^  <  0  ;  if  a^  ■  0  then  neither  u^  nor  appears 

in  . 
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An  assignment  for  the  variables  u^,...,un  is  a  function  o  :  {uj,..., 
uq}  -*■  {F,  T}  .  A  clause  Is  satisfied  by  a  If  there  exists  j  such 

that  Uj  appears  In  and  a(u^)  -  T  or  Uj  appears  In  and  a(uj)  * 

F  .  The  formula  B  Is  satisfied  by  a  If  each  clause  In  B  Is  satis¬ 

fied  by  a  .  The  satisfiability  problem  Is  that  of  recognizing  formulas  B 
that  are  satlsflable.  NP-completeness  of  this  problem  was  established  by 
Cook  [4],  and  that  was  the  fundamental  result  from  which  the  theory  of  NP- 
completeness  was  developed  by  Karp  [14]  and  others  (see  Garey  and  Johnson  [10] 
for  an  extensive  survey) .  It  seems  very  unlikely  that  the  general  satisfia¬ 
bility  problem  admits  a  polynomlally  bounded  algorithm. 

A  weak  assignment  for  the  variables  Uj,...,un  is  a  function  a  :  {uj, 
...,un>  (F,  Z,  T}  (where  Z  may  be  regarded  as  standing  for  "zero")  such 

that  for  at  least  one  j  ,  a(u^)  #  Z  .  A  clause  Is  weakly  satisfied 

by  a  If  Is  satisfied  by  a  or  a  misses  in  the  sense  that  a(Uj)  “ 

Z  for  all  j  such  that  Uj  or  Uj  appears  In  .  The  formula  B  is 
weakly  satisfied  by  a  If  each  Is  weakly  satisfied.  The  weak  satisfia¬ 

bility  problem  Is  that  of  recognizing  formulas  B  that  are  weakly  satlsflable. 
Plainly  a  matrix  A  Is  weakly  satlsflable  (that  is,  A  Is  a  W-matrix)  if 
and  only  If  the  corresponding  Boolean  formula  B  Is  weakly  satlsflable. 

Three  examples  serve  to  illustrate  the  above  notions. 

3x2  matrix  Boolean  Formula  Properties 

♦  0  (Uj)  A  (Uj  vu2)  a  (Uj  V  u2)  Satisfied  by  o(uJ)  ■  T  , 

♦  -  a(Uj)  »  F  .  Weakly  satisfied 

♦  -  (but  not  satisfied)  by  a(u^)  ■ 

Z  »  a(u2)  ■  F  . 


i 


1 


3  x  2  matrix 


Boolean  Formula 


(up  a  (u2)  a  (up 


Properties 


Not  satisfiable.  Weakly  satis¬ 
fied  by  a(up  «  Z  ,  a(up  *  T 


(u1  v  u2)  a  (up  v  (G2) 


Not  weakly  satisfiable. 


We  now  show  how  to  reduce  the  satisfiability  problem  to  the  weak  satis¬ 
fiability  problem. 

THEOREM  A  Suppose  that  . are  clauses  in  the  n  variables  u0,..., 

u  .  and  D,,...,D  are  clauses  in  the  2n  variables  un . u  v_,..., 

n-  X  x  m  u  n- x  u 

vn  j  ,  where  is  obtained  from  by  replacing  each  occurrence  of 

with  Vj  .  Form  An  additional  clauses  Ej ,  Fj,  Gj  and  as  follows , 
where  the  subscript  j  ranges  from  0  to  m  -  1  and  j  +  1  is  reduced 
modulo  m  : 

Ei  iB  uj  v  vj  •  fj  iB  zs  v  • 

Gj  ie  v  Vi v  V»  •  "j  u  ’3 v  “3+1 v  Vi  • 

For  the  formula  B  •  Cj  a  ...  a  c  and  the  formula 

B'  -  D,  A  ...  A  D  A  E.  A  ...  A  E  A  F,  A  ...  A  F_  A 
1  ml  n  1  n 

G.  A  ...  aG  aH.  a  ...  AH  ■ 

1  n  1  n 

the  following  three  conditions  are  equivalent:  B  is  satisfiable ;  B'  is 
satisfiable i  B*  is  weakly  satisfiable. 

Proof.  If  B  is  satisfied  by  an  assignment  o  :  {u^,...^}  -*■  (F,  T) 
then  B*  is  satisfied  by  the  assignment  o'  :  (uj,...,u  ,  vi»“*»vn* 

(F,  T}  ,  where  a’(up  -  o(Uj)  and  a'(vp  -  “*(up  .  Conversely,  if  a' 
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ls  any  assignment  that  satisfies  B'  then  the  restriction  of  a'  to 
{u, . u  }  satisfies  B  because  from  the  satisfaction  of  E.  and  F, 

1  n  J  j 

if  follows  that  a*  assigns  complementary  values  to  u^  and  . 

To  complete  the  proof  we  show  that  if  8  :  {uj,...,un,  Vj . vr}  -*• 

{F,  Z,  T}  is  a  weak  assignment  that  weakly  satisfies  B'  then  Z  is 
missing  from  the  range  of  8  ,  whence  8  is  in  fact  an  assignment  that 
satisfies  B'  . 

By  the  definition  of  weak  assignment,  there  exists  k  such  that 
SCuj^)  4  Z  or  8(vk)  4  Z.  For  any  such  k  ,  if  BCv^)  *  T  then  8(vk>  *  F 
by  Fk  ,  while  if  8^)  •  F  then  6(vk>  *  T  by  Efc  .  Similarly,  if 
B(vk)  -  T  then  P(uk>  *  F  by  Ffc  ,  while  if  6(vk>  *  F  then  B(uk>  «  T 

by  V 


It  remains  only  to  show,  by  Induction  on  J  ,  that  if  BC^)  e  {T,  F} 
then  BCu^j)  e  {T,  F}  for  all  j  .  But  if  e(uk+^)  ■  T  then  by 
one  of  811,1  oust  be  assigned  T  ,  whence  the  other  is 

assigned  F  by  the  observation  of  the  preceding  paragraph.  And  if  8 ( uk+j ) 
F  then  B(vfe+j)  -  T  ,  whence  by  one  of  ufc+j+1  and  vfc+j+1  must  be 

assigned  T  ;  hence  the  other  is  assigned  F  .  That  completes  the  Induction 
and  the  proof.  □ 

For  each  k  >  0  ,  let  the  decision  problem  WSATk  be  as  follows: 

Instance:  A  Boolean  formula  B  in  conjunctive  normal  form,  involving 

a  total  of  n  propositional  variables  and  consisting  of  the  conjunction  of 
1/k 

n  +  Ln  J  disjunctive  clauses; 

Question:  Is  B  weakly  satisfiable? 
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THEOREM  5  For  each  k  >  0  the  problem  WSAT^  is  NP -complete. 

Proof.  Plainly  WSAT^  belongs  to  NP  ,  so  it  suffices  to  describe 
a  polynomial  reduction  of  SAT  (the  usual  satisfiability  problem)  to 
WSAT^  .  Consider  an  instance  of  SAT  that  concerns  a  Boolean  formula 
B  in  conjunctive  normal  form,  involving  a  total  of  n  propositional 
variables  and  consisting  of  the  conjunction  of  m  disjunctive  clauses. 

Let  the  variables  be  uo»***»un_i  »  int*°duce  additional  variables  v  , ..., 

If 

,  and  let  the  formula  B'  be  as  in  Theorem  4.  Let  r  *  (m  +  2n  +  1)  , 

introduce  additional  variables  Wg,...,vr_j  ,  and  let  B*  be  a  not-weakly- 
satisfiable  formula  consisting  of  the  conjunction  of  r  +  1  disjunctive 
clauses  in  these  variables  (easily  constructed  with  the  aid  of  Theorem  1) . 
Finally,  let  B"  denote  the  conjunction  of  B'  and  B*  .  Then  B"  con¬ 
sists  of  m  +  4n  +  r  +  1  clauses  and  involves  a  total  of  2n  +  r  variables. 
Since  the  inequality 

m+4n+r+l<  (2n  +  r)  +  (2n  +  r)*^ 
is  equivalent  to  the  valid  inequality 
(m  +  2n  +  1)^  <  2n  +  r  , 

some  of  the  clauses  of  B"  can  be  repeated  if  necessary  to  produce  an 

instance  Bm  of  WSAT,  .  From  Theorem  4  and  the  fact  that  B*  is  not 

k 

weakly  satisfiable,  it  follows  that  B  is  satisflable  if  and  only  if  Bm 
is  weakly  satisfiable.  □ 

As  can  be  seen  from  the  references  in  Garey  and  Johnson  [10],  there 
are  several  other  relatives  of  the  satisfiability  problem  for  which  NP- 
completeness  has  been  established.  See  especially  Schaefer  [30]. 
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4.  The  Recognition  of  S-svstems 


„m 


As  the  term  Is  used  here,  a  cone  Is  a  nonempty  subset  C  of  R 
such  that 

0  k  C  -  C  +  C  -  30,»[C  . 

Thus  C  omits  the  origin  and  is  closed  under  vector  addition  and  multi¬ 
plication  by  positive  scalars;  in  particular,  C  is  convex.  A  cone  is 
■polyhedral  if  it  is  defined  by  a  finite  system  of  linear  inequalities. 

An  S -system  is  a  sequence  (C^, . . . ,C  t  of  cones  in  Rm  such  that  for 

each  m  x  (mfl)  matrix  B  whose  column  belongs  to  (1  <  k  £  m  +  1) 

m| 

and  for  each  x  e  R  with  Bx  *  0  ,  all  coordinates  of  x  are  of  the  same 
sign.  Thus  an  m  x  (m+-l)  matrix  with  columns  vi>***>vm+i  is  311  S-matrix 
if  and  only  if  the  sequence  (Q(Vj) , . . .  j))  is  an  S-system. 

Despite  the  NP-completeness  result  of  Section  3  and  the  close  relation¬ 
ship  between  S-matrices  and  weak  satisfiability,  it  can  be  determined  in 
polynomial  time  whether  a  given  sequence  of  polyhedral  cones  is  an  S-system. 
The  algorithm  is  based  on  the  following  characterization. 

THEOREM  6  If  C^. . . . ,Cm>  ^  are  cones  in  Rm  and  e  for  each  k  , 
then  the  following  four  conditions  are  equivalent: 

(i)  (Cj,. ie  an  S-system i 

(ii)  for  each  choice  of  Cj  c  »  the  c^'s  are  the 

vertices  of  an  m-s implex  whose  interior  contains  the  origin ; 

(ill)  for  each  nonzero  z  e  Rm  there  exists  k  such  that  cCz  >  0 
for  all  c  e  ; 

(iv)  for  each  choice  of  distinct  r  ,  s  e  {1,. . . ,®fl)  the  linear  hull 
of  the  set  {p^  :  r  +  k  «}  ie  a  hyperplane  that  strictly  sepa¬ 
ratee  Cr  from  Cg  . 
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Proof.  Conditions  (i)-(iii)  are  merely  the  extensions,  to  general 
cones,  of  the  corresponding  conditions  of  Theorem  2.  The  reasoning  of 
Theorem  2  applies  here  as  well.  Our  algorithm  is  based  on  the  equivalence 
between  (ii)  and  (iv),  which  will  now  be  established. 


If  (ii)  holds  and  c^  e  for  each  k  then  the  points  c^  are 
affinely  independent  and  the  affine  hull  of  any  m  of  them  is  a  hyperplane 
in  Rm  that  misses  the  origin.  This  implies  that  each  m  of  the  c^'s 
form  a  linearly  independent  set,  whence  the  linear  hull  of  any  m  -  1  of 
them  is  a  hyperplane  through  the  origin.  The  remaining  two  c^’s  are  on 
opposite  sides  of  that  hyperplane,  for  the  origin  is  interior  to  the  convex 
hull  of  the  c^'s  .  It  follows  that  (ii)  implies  (iv). 

Now  suppose  that  (iv)  holds,  let  T  denote  the  set  of  all  (nri-1) -tuples 
(Cj, . . . ,  c^)  such  that  c^  e  for  all  k  ,  and  let  r*  denote  the  set 
of  all  (m  +  l)-tuples  in  T  such  that  for  each  choice  of  distinct  r  , 


s  e  {l,...,m  +  1}  the  linear  hull  of  the  set  (c^ir^k^s)  is  a  hyper¬ 
plane  that  strictly  separates  C^.  from  Cg  .  Then  (p^ ,  • . .  tP^ j)  t  r*  by 
(iv),  and  we  claim  that  in  fact  T*  ■  T  .  To  prove  this  it  suffices  to  show 


that  if  (c^, . . .  ,c  j  j)  er*,ls.£<m+l,  and  is  the  set  of  all 

d  e  such  that  T*  Includes  the  (m  +  l)-tuple  obtained  from  (c^,..., 
c|m)  p  by  replacing  c ^  with  d  ,  then  . 

For  each  e  e  and  each  r  e  (l,...,m  +  1}  ~  {£}  let 
Xr(e)  -  (e)  u  {(^  :  k  k  U,  r}}, 
and  for  each  s  e  {l,...,m  +  1}  ~  {r,  s)  let 


Xrg(e)  »  {e}  u  {c^  :  k  4  U,  r,  s}}  . 

The  set  Xr(e)  is  linearly  Independent  because  the  relevant  c^'s  are  linearly 
Independent  and  their  linear  hull  misses  e  .  Hence  the  set  Xrg(e)  is  linearly 
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independent  and  Its  linear  hull  is  a  hyperplane  H  (e)  through  the 

rs 

origin.  To  show  e  e  it  remains  to  show  that  Hrg(e)  strictly 

separates  Cr  from  Cg  .  Since  c£  e  »  such  separation  does  occur 

when  e  e  ]0,®Cc£  .  Hence  if  the  desired  separation  fails  to  occur 

for  the  e  in  question,  there  exists  d  £  ]  c^,e  [  ~  ]0,®[c^  such  that 

HfS(d)  Intersects  C  or  Cg  .  We  suppose  without  loss  of  generality 

that  Hrg(d)  Intersects  Cr  .  Then  there  exists  b  e  Cr  and  there 

exist  scalars  A,  such  that 
k 

b  “  XlA  +  hkU,T,s)  * 

Also,  since  the  linear  hull  of  {c^  :  r  4  k  4  s}  is  a  hyperplane  that 
strictly  separates  Cr  from  Cg  ,  there  are  scalars  such  that 

b  -  vtct  *  IW{riS)  ^  . 

Plainly  A^  <  0  and  <  0  .  Since  the  set  Xr(b)  is  linearly  inde¬ 
pendent  it  follows  from  the  two  equations  that  d  e  ]t3,“[c£  ,  a  contra¬ 
diction  showing  that  .  It  follows,  then,  that  T*  *  r  . 

Now  consider  an  arbitrary  (m  +  l)-tuple  (Cj , . . .  ,0^)  c  T  ■  T* 
and  an  arbitrary  choice  of  scalars  y ^  such  that  *  0  •  If 

there  exist  r  and  s  such  that  y  s  0  £  y  then  all  are  0  , 

as  follows  from  the  equation 

V.  *  (^*r  -  IrfkfcVk 

in  conjunction  with  the  fact  that  the  linear  hull  of  (c^  :  r  O  ^  is 
a  hyperplane  that  strictly  separates  cr  from  cg  .  This  condition  on  the 
Yjj's  implies  the  c^’s  ere  affinely  Independent  and  are  in  fact  the  ver¬ 
tices  of  an  m-simplex  whose  interior  contains  the  origin.  Hence  (iv)  implies 
(ii).  0 
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The  following  theorem  applies  to  arbitrary  cones  but  may  be  difficult 
to  Implement  if  the  cones  are  not  polyhedral. 

THEOREM  7  The  following  procedure  decides  correctly  whether  a  given  sequence 
(Cj, . . .  j)  of  cones  in  Rm  is  an  S-system: 

(1)  Select  points  p1  e  Ci»***»PnH-l  £  Cnrt-1  ' 

(2)  Decide  whether  there  exists  a  nonzero  z  e  Rm  such  that  p*z  >  0 

for  all  k  .  If  so j  stop ;  (C^ . CnH-p  ^ B  noi  an  S-system. 

(3)  For  each  r,s  with  l£r<s£m+l  ,  form  the  linear  hull  H 

ITS 

of  (p,  :  r  i*  k  ^  s}  .  Decide  whether  H  intersects  C  and 
K  ITS  r 

whether  Hrg  intersects  Cg  .  If  a  nonempty  intersection  is 
foundt  stop;  (Cj, . . .  j)  is  not  an  S-system  . 

(4)  If  Cr  n  Hrg  *  0  «  Cg  n  Hrg  for  all  r  <  s  then  (C^, . . .  .C^ x) 
is  an  S-system  . 


Proof.  If  there  exists  a  z  as  in  (2)  then  {p^ . p^ lies  in  a 

closed  halfspace  whose  bounding  hyperplane  passes  through  the  origin.  Hence 


the  origin  is  not  Interior  to  the  convex  hull  of  the  p^’s  and  (C^, . . .  .C^) 
is  not  an  S-system.  If  no  such  z  exists  then  each  m  of  the  points  p^,..., 
Pm) i  are  linearly  independent  and  each  Hrg  is  a  hyperplane.  If  Hrg  inter¬ 
sects  C  or  C  ,  we  see  as  in  the  proof  of  Theorem  6  that  (C . C.  . .)  is 

not  an  S-system.  But  if,  in  step  (3),  the  hyperplane  Hfg  misses  both  Cr 
and  C  then  It  strictly  separates  them;  for  otherwise  C  and  C  lie  to- 

gather  In  a  closed  half space  bounded  by  Hrg  ,  and  a  z  would  have  been  found 
(and  the  computation  halted)  In  step  (2) .  Thus  if  the  confutation  progresses 


past  step  (3),  it  follows  from  Theorem  6  that  (C^,. . .  is  an  S-system.  □ 
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When  the  cones  are  polyhedral,  the  selection  of  the  p^'s 

In  step  (1)  and  the  consistency  tests  in  steps  (2)  and  (3)  can  be  car¬ 
ried  out  by  means  of  various  linear  programming  methods.  In  particular. 

Phase  I  of  Dantzig's  simplex  method  [5]  is  applicable,  as  in  the  recent 
Shor-Khachian  method  [15,  31].  In  view  of  numerous  reports  of  computa¬ 
tional  experience  with  these  two  methods  [3,  5,  6,  25,  32],  we  may  say 
that  each  leads  to  a  "good"  algorithm  for  recognizing  polyhedral  S-systems. 
The  recognition  algorithm  based  on  the  simplex  method  is  good  in  the  sense 
that  it  is  usually  very  efficient  in  practice,  even  though  its  worst-case 
behavior  is  not  polynomially  bounded  [16].  The  recognition  algorithm  based 
on  the  Shor-Khachian  method  is  good  in  the  sense  that  it  is  polynomially 
bounded  (see  Theorem  8),  even  though  it  is  usually  very  inefficient  in 
practice.  It  would  be  very  nice,  for  the  recognition  of  polyhedral  S-systems 
as  well  as  for  solving  linear  programming  problems,  to  have  an  algorithm  that 
is  good  in  both  senses! 

THEOREM  8  By  using  the  procedure  of  Theorem  7  in  conjunction  with  the 
Russian  method  for  linear  inequalities ,  the  problem  of  recognizing  polyhedral 
S-systems  can  be  solved  in  polynomial  time  and  space  on  a  deterministic  Tur¬ 
ing  machine.  Specifically ,  when  each  of  the  cones  Cj . C^j  in  Rm  is 

given  by  a  system  of  linear  inequalities  with  integer  coefficients  and  E 
is  the  length  of  the  binary  encoding  of  the  input  datat  the  algorithm  decides 
in  time  0(m5E)  whether  (CJt . . ie  an  S -system. 

Proof.  It  follows  from  Khachlan's  analysis  [15]  of  the  Shor-Khachian 
algorithm  [15,  31]  that  a  polynomially  bounded  algorithm  for  recognizing 
polyhedral  S-systems  Is  obtained  If  one  implements  the  procedure  of  Theorem  7 
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by  suitable  use  of  the  SK-algorithm  to  find  the  points  pk  and  to  make 
the  required  consistency  tests.  Khachlan's  paper  [15]  is  an  extended 
abstract,  without  proofs,  but  proofs  and  various  Improvements  have 
appeared  in  technical  reports  Issued  subsequently  by  many  operations 
researchers,  computer  scientists  and  mathematicians.  We  mention  Gdcs  and 
Lovdsz  [8],  Aspvall  and  Stone  [1],  and  especially  Padberg  and  Rao  [26,  27] 
who  have  an  improvement  of  the  SK-algorithm  that  they  call  the  Russian 
method  for  linear  inequalities.  Our  discussion  is  based  on  that  of  Padberg 
and  Rao. 

In  order  to  handle  both  strong  and  weak  linear  inequalities,  let  us 
assume  that  for  each  k  there  are  finite  sets  and  of  nonzero 

lattice  points  in  Rm  such  that 

-  {x  e  Rm  :  yCx  >  0  for  all  y  e  Y^  and  ztx  2  0  for 

all  z  c  Zfc}  . 


The  set  may  be  empty,  but  we  may  assume  Y^  is  nonempty  for  otherwise 

0  c  and  (Cj, . .  .  .C^j)  is  not  an  S-system.  Equality  constraints  w*x  *  0 
are  handled,  as  usual,  by  including  both  w  and  -w  in  the  set  .  In 
particular,'  if  is  of  the  form  Q(V^)  ,  {u^,...,um)  is  the  standard  basis 

for  R01  ,  and  L  of  V^’s  coordinates  are  0  ,  then  Y^  includes  u^  or 
-u^  for  m  -  l  values  of  1  and  includes  both  u^  and  -u^  for  the 
remaining  t  values  of  1  . 

k  “k 

For  each  k  let  n^  ■  lY^I  +  IZ^I  ,  let  b  be  the  point  of  R 

whose  first  lY^I  coordinates  are  1  and  last  |Z^|  coordinates  are  0  , 
k  k 

let  A  ■  (*£j)  be  the  n^  *  m  matrix  whose  first  lY^I  rows  are  the 


vectors  y*  for  y  c  and  last  IZ^I  rows  are  the  vectors  z* 
for  z  €  ,  and  let 

-  {x  e  Rm  :  Akx  2:  bk}  c  . 

Note  that  a  linear  subspace  of  KSi‘  Intersects  if  and  only  if  it 

Intersects  .  Thus  It  is  permissible.  In  following  the  procedure 
described  in  Theorem  7,  to  work  with  the  Dk's  rather  than  the  C^'s  . 

If  we  adopt  the  special  convention  that  log^O  -  0  ,  the  length  of 
the  binary  encoding  of  the  description  of  (or  of  )  may  be  taken 

as 

Ek  -  1  +  Llog2nU  +  1^  , 

where 

1^  -  2  +  llog2|YkU  +  Llog^J  +  200^  +  Jj.j  Llog2 |a^j  I J  . 

Here  I |  denotes  cardinality  and  ta^jl  denotes  absolute  value.  We 
assume  that  for  each  integer  x  ,  1  +  log2Lxj  bits  are  used  to  represent 
| x I  and  an  additional  bit  represents  the  sign  of  x  if  the  problem's 
formulation  does  not  automatically  require  x  >  0  (as  it  does  when  x  is 
m  or  nk  ) .  The  length  of  the  encoding  of  the  description  of  the  system 
is  defined  as 

E  -  1  +  Llog2mJ  +  1^  . 

It  is  easy  to  quibble  in  minor  ways  with  these  definitions  of  the  Ek's 
and  of  E  ,  but  such  quibbles  do  not  affect  the  validity  of  our  theorem. 
Let 

Ffc  -  1  +  llog2mJ  +  ^  , 

where  Ak  is  the  maximum  of  the  absolute  values  of  the  determinants  of 


k  k 

the  square  submatrices  of  the  matrix  [A  ,  b  ]  ,  and  let 

.  F. 

D£  *  {x  £  :  |x^|  s  2^  2  for  all  1  }  . 

It  follows  from  the  discussion  of  Padberg  and  Rao  [26,  27]  that  , 

nonemptiness  of  Implies  nonemptiness  of  ,  and  the  Russian  method 

3 

yields  a  point  pfc  e  In  time  (and  space)  0(m  F^)  .  Thus  the  complexity 
of  step  (1)  in  the  procedure  of  Theorem  7  is 

0<m3l£jFk)  5  0(m3E)  . 

With  pfc  e  D£  for  1  £  k  £  m  +  1  ,  the  consistency  problem  of  step  (2) 

admits  a  binary  encoding  of  length  0(E)  and  hence  step  (2)  Is  handled  in 

time  OOn^E)  by  the  Russian  method. 

To  form  the  linear  hull  H  of  (p.  :  r  +  k  i  s}  ,  as  required  by 

rs  k 

step  (3),  let  M  be  the  (m  -  1)  x  m  matrix  whose  columns  are  the  p.'s 
rs  k 

for  r  i*  k  j*  8  ;  then  use  the  Russian  method  (or  elementary  column  operations) 

to  find  a  nonzero  q  c  Rm  such  that  q^  M_  *  0  and  all  coordinates  of 

rs  rs  rs 

3 

q  are  less  than  E  In  absolute  value.  This  can  be  done  In  time  0(m  E)  , 

>*rs 

and  then 


-  <*  «  *"  :  4s* 


rs 


0}  . 


Testing  whether  H  Intersects  C  or  C  Is  also  of  complexity  0(m3E) 
rs  r  s 

If  the  Russian  method  is  used.  Since  there  are  m(m  +  l)/2  pairs  (r,  s) 
to  contend  with,  the  overall  complexity  of  step  (3)  and  of  the  entire  algo¬ 


rithm  is  0(n>  E)  .  □ 
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5.  Open  Problems 

If  each  (n  +  1)  *  n  NW-matrix  can  be  generated  from  smaller  ones  by 
the  method  of  Theorem  2,  then  each  S-matrix  has  at  least  one  row  in  which 
there  are  at  most  two  nonzero  entries.  Gorman  has  conjectured  that  this  is 
so  and  can  be  used  as  the  basis  of  a  fast  algorithm  for  recognizing  S-matrices. 

Let  *(m,  n)  denote  the  number  of  equivalence  classes  of  m  x  n  NW- 
matrices.  What  can  be  said  about  the  asymptotic  behavior  of  the  function  *  , 
or  about  its  values  for  small  (m,  n)  ?  There  is  special  interest  in 
*(n  +  1,  n)  ,  the  number  of  equivalence  classes  of  n  *  (n  +  1)  S-matrices. 

We  have  6een  that  #(2,1)  *  1,  *(3,2)  «  2,  *(4,3)  *  10  .  What  is  *(5,4)  ? 

Consider  the  problem  of  recognizing  whether  a  given  m  *  n  matrix  is 

an  NW-matrix.  As  we  have  seen,  this  problem  can  be  solved  in  polynomial  time 

1/k 

when  m  ■  n  +  1  but  is  NP-complete  when  m  -  n  +  In  \J  for  an  arbitrary 
fixed  k  >  0  .  What  happens  when  m  ■  n  +  k  for  a  fixed  k  >  1  ? 

For  each  n  let  h(n)  denote  the  smallest  integer  r  that  has  the 
following  property:  whenever  a  Boolean  formula  is  in  conjunctive  normal 
form.  Involves  a  total  of  n  propositional  variables,  and  is  not  weakly 
satlsfiable,  then  some  subformula,  consisting  of  the  conjunction  of  h(n) 
or  fewer  of  the  clauses  of  the  original  formula,  is  not  weakly  satlsfiable. 

It  is  easy  to  see  that  2n  £  h(n)  £  3n  -  1  .  Plainly  h(l)  •  2  and  it  can 
be  verified  that  h(2)  ■  4  .  Is  h(n)  always  equal  to  2n  ? 

For  each  x  e  Rn  ~  {0}  ,  let  U(x)  denote  the  set  of  all  unit  vectors 
u  e  Rn  such  that  u  weakly  satisfies  c  .  Since  U(x)  depends  only  Q(x)  , 
there  ere  only  3n  -  1  sets  of  the  form  U(x)  for  x  e  Rn  -  {0}  .  For 
example,  when  n  ■  2  there  ere  four  sets  like  each  of  the  samples  shown  in 
Fig.  9. 
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2 

Fig.  9:  The  two  types  of  sets  U(x)  for  x  e  R  ~  {0} 

Let  us  use  the  term  U  - set  to  denote  a  set  of  the  form  U(x)  for 

n 

x  €  Rn  ~  {0}  .  Then  h(n)  is  the  smallest  k  such  that  whenever  each 
k  members  of  a  family  of  l^-sets  have  a  common  point,  the  entire 
family  has  nonempty  Intersection.  Thus  h(n)  is  a  Helly  number  in  the 
sense  of  Danzer,  Grunbaum  and  Klee  [73.  It  seems  possible  that  methods 
similar  to  those  of  Grunbaum  and  Motzkin  [13]  would  be  useful  in  deter¬ 
mining  h(n)  . 

For  each  cone  K  in  Rn  ,  let  W(K)  denote  the  set  of  all 
w  €  Rn  ~  {0}  such  that  w*k  2:  0  for  each  k  c  K  .  For  each  family 
£  of  cones,  let  g(£)  -  (W(K)  s  K  e  &}  and  let  H(jK)  denote  the  Helly 
number  of  W(K)  .  Then  h(n)  is  H({Q(x)  :  x  e  Rn  ~  {0})  .  What  can  be 
said  about  H(K)  for  other  interesting  choices  of  K  ? 
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